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^ ■ Abstract 

^ , Random matrix theory (RMT) provides a common mathematical formulation 

(~| , of distinct physical questions in three different areas: quantum chaos, the 1-d 



integrable model with the interaction (the Calogero-Sutherland-Moser 
system), and 2-d quantum gravity. We review the connection of RMT with 
these areas. We also discuss the method of loop equations for determining cor- 
relation functions in RMT, and smoothed global eigenvalue correlators in the 
2-matrix model for gaussian orthogonal, unitary and symplectic ensembles. 

1. Introduction 

In the last few years it has been realized that matrix models form a bridge between two 
apparently distinct subjects: the Calogero-Sutherland-Moser integrable model of particles 
in 1-d interacting via the potential (CSM model), and quantum chaos. Correlations of 
the eigenvalue density in quantum chaotic systems and spacetime dependent particle density 



1 



correlators in the CSM model both map onto the same eigenvalue correlator in a random 
matrix model. Furthermore very closely related correlation functions in the matrix model 
appear in the context of QCD, string theory and 2-d quantum gravity - the correlation 
functions of 'loop operators', which have a geometric interpretation in terms of surfaces. 
Thus different physical questions in these three vastly different areas reduce to the same 
mathematical question formulated in terms of matrix models. In this review we discuss these 
questions and their common mathematical formulation. We also review the method of loop 
equations, a technique for calculating these correlators that originated in the QCD/string 
theory context, and its application to obtain correlation functions in the 2-matrix model 
that are relevant for quantum chaos and the CSM model. 

In section 2 an introduction to RMT is given. We collect some recent results for cor- 
relation functions in 1- and 2-matrix models that are relevant to the sequel. Sections 3 
and 4 describe the apphcation of RMT in the field of quantum chaos and for obtaining the 
static correlators of the CSM model, respectively. Section 5 discusses a deeper relationship 
between quantum chaos and the CSM model, that is, the equivalence of parametric correla- 
tions in quantum chaos and time dependent correlations in the CSM model via a mapping 
onto a 2-matrix model. Two-dimensional quantum gravity and its relationship with matrix 
models is described in section 6. Section 7 discusses loop equations for determining corre- 
lation functions in matrix models and derives some of the results for eigenvalue correlation 
functions mentioned in section 2. Section 8 contains a brief summary. 

2. Random Matrix Models 

Consider the partition function 



where ^4 is an A'^ x hermitian matrix, dA is the standard U{N) invariant measure for 
hermitian matrices, and S{A) is a. U (N) invariant action of the form 

OO -I 

SiA) = NTrViA), V{A) = ^ -g^A\ (2) 

n=l ^ 




(1) 
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This partition function defines the normahzation of a probabihty distribution for the matrix 
A, namely P{A) = (1/Z)e'~^^^\ in which the expectation value of any function f{A) is 
given by (/) = / dA P{A) f{A). Typical examples of functions whose expectation values 
and correlation functions will be of interest are the Green's function or 'loop operator' 
W{z) = {1/N)Tt{z — A)^^, and the eigenvalue density operator p{x) = (l/A^)Tr 6{x — A). 
One is mostly interested in this matrix model in the large N limit. 

The U{N) symmetry of the above probability distribution (namely symmetry under the 
transformation A — > UAU'^, with U any N x N unitary matrix) can be used to project it 
down to a joint probability distribution of the eigenvalues Xi, Xn of A. This probability 
distribution is given by 

P(Ai, A^) = iVo n |A. - X,f e'^ v^ca.)^ (3) 

where (3 = 2. The projection is achieved by using the symmetry to diagonalize A : 
A = UDU^, D = diag(Ai, Aiv). The 'angular variables' U decouple from the eigen- 
value integration because of the U{N) symmetry and their integration provides merely an 
overall factor absorbed in the normalization constant Nq. The Van der Monde determinant 
factor Yli^j \Xi — Xjf is simply the jacobian of the change of variables from A to D and U. 

As a consequence of this projection the expectation value of any 'observable' (any U (N) 
invariant function f{A); in particular W{z) and p{x) above) is given by 

(/) = iVo 1 dX,...dXN e-^(^^--'"-) fiD), (4) 

with 

N 

E{X,, Xn) = NJ2 V{X,) -PY.hg |A, - A,|. (5) 

1=1 i<j 

This provides a very useful physical picture of the matrix model: the eigenvalue A, can be 
thought of as the position coordinate of the ith particle in a 1-dimensional gas of particles 
(the Dyson gas |T]]) in which all particles experience the same external potential NV{x), and 
repel each other logarithmically, giving rise to a (static) energy function E. (/) is merely 
the Boltzmann weighted average of the function f{D) of the position variables. 

3 



If A belonged to the ensemble of N x N real symmetric matrices rather than hermitian 
matrices (that is to say, for real symmetric A, P{A) is the same as that given above, but 
zero otherwise), equations (0-|^) remain unchanged, except that now (3=1. Similarly, 
if A belonged to the ensemble of N x N real self-dual quaternions, then j3 = 4. The 
three ensembles corresponding to (3 = 1,2,4 are referred to as the orthogonal, unitary 
and symplectic ensembles respectively, corresponding to the symmetry group in each case 
For details, see, e.g., [§]. /?, which determines the strength of the eigenvalue repulsion, 
simply measures the unlikelyhood of finding a matrix in the ensemble with a pair of equal 
eigenvalues. 

The probability distribution (|) depends upon the parameters (?„ of the potential V. The 
simplest cases are the gaussian ensembles (where g2 = fi > 0, gi = gn = for n > 3) for 
which the one-point function (or eigenvalue density) p{x) = (p(x)) obeys the semicircle law 



with = 2j3/fi. In general, when V{x) has a single minimum, p{x) is nonvanishing on a 



single segment [a,b] of the real line and has the form p{x) = P{x)\J{x — a){b — x), where 
P{x) is a polynomial in x whose degree depends upon the degree of V. Both a and b as well 
as the coefficients of P{x) depend upon the gn, thereby making p{x) nonuniversal. 

However, it turns out that higher correlation functions of the eigenvalue density and 
the loop operator exhibit a remarkable universality (independence of the choice of gn) in 
certain domains. For example, the exact two point density-density correlator {p{x)p{y)) for 
an arbitrary even polynomial potential V for the unitary ensemble {jS = 2), computed in 
0, though not universal in general, is universal in two physically interesting domains. The 
first domain is the 'local' region where the interval [x,y] contains a finite number {0{N^)) 
of eigenvalues and is located a finite distance from the endpoints. Then, the result, for 
arbitrary V, is 



i 



p{x) = — -Va^ — x"^, \x\ < a, and p{x) = 0, |x| > a, 



(6) 




{p{x)p{y)) 
p{x)p{y) 



1 - 




5 = nN{x — y)p{ 



X + y 
2 



(7) 
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which is universal when expressed in terms of variables that are rescaled by the factor 
P(^2^) ^hat determines the local mean level spacing. This is a 'fine grained' correlator that 
oscillates over distances of the order of mean eigenvalue spacing, and has a nonsingular limit 
a.s X — y 0. was derived originally for the circular unitary ensemble (CUE: ensemble 
of unitary matrices, whose eigenvalues lie on a circle) in ^ and for the gaussian unitary 
ensemble (GUE: ensemble of hermitian matrices with a quadratic V^) in 0. Its universality 
for the unitary ensemble (/3 = 2) was extended to all V^s that correspond to the classical 
orthogonal polynomials in p|] and to all even polynomial V^'s in Expressions analogous 
to d^) exist for the other two ensembles (/5 = 1,4) The universality of all these local 



fine grained correlators with respect to arbitrary potentials has been shown in |jTO[- For 



other extensions of the universality of (0) and the higher point local correlators see [pT 



While the universality of local correlators has long been expected on empirical grounds 
(see next section), it is somewhat more of a surprise that even certain 'global' correlators are 
universal. This is the domain where [x, y] contains a large number of eigenvalues including 
upto 0{N), but the correlator is 'smoothed' by independently averaging over x and y over 
intervals much greater than the mean eigenvalue spacing. In this domain the result for the 
connected two point function in the large limit is 

This is universal in that it depends on the Qn only through the endpoint of the spectrum 
a. (For even potentials the endpoints are ±a. For non-even potentials, the result is still 
(I), but now a represents half the width of the support, and x,y on the r.h.s. of (||) are 
shifted by the midpoint of the support.) This is a 'global' correlator in the sense that x and 
y could be anywhere in the interval [—a, a], close to the endpoints or far away (provided 
only that x — y contains a large enough number of eigenvalues for 'smoothing'). (M) was 



derived originally for the gaussian ensembles in ref. |12|. Its universality with respect to 



arbitrary potentials for /5 = 2 follows from the results of ref. |TB| obtained for the two-point 
function of loop operators. In |^ (^ was shown to be the smoothed version of the exact 



two-point density- density correlator for arbitrary even polynomial potentials in the unitary 
ensemble. (^) has been proven to be universal for arbitrary potentials V which give rise to 



an eigenvalue density with support in a single segment for all three ensembles [0]. For these 
potentials the three and higher point smoothed correlators vanish to leading order in 
P JT^JT3|] and are in general nonuniversal in higher orders |TB . 



There is yet another type of universality of smoothed correlators at multicritical points 
1 171 that arises in the double scaling limit [jl8l studied in the context of 2-d gravity and 



noncritical string theory [|T^]. Here only the region close to the endpoint of the support 
of p{x) matters, and correlations depend on the Qn only through the number of zeroes of 
the polynomial P{x) that coincide with the endpoint a (the level of multicriticality) ||20|| . 
Further, there is evidence of a universality of finegrained correlators at the spectrum edge 



even in the absence of a double scaling limit [21 



The 2-matrix model is defined by the partition function 

dA (IB e"^ (9) 



where A and B are N x N matrices drawn from any of the above mentioned ensembles, and 
we will restrict the discussion to S of the type S = NTi (Vii^A) + V2{B) — cAB). When A 
and B are drawn from different ensembles, e.g., A is drawn from the gaussian orthogonal 
or symplectic ensemble (/3 = 1 or 4) and B from the gaussian unitary ensemble {jS = 2), 
the eigenvalue correlations (relevant for parametrizing time reversal symmetry breaking in 



chaotic systems) have been computed in Here we will be primarily interested in the 



case where both A and B are drawn from the same ensemble. 

The smoothed, global expression for the 2-point function pAB{x,y) = {pA{x)pB{y))c = 
{pA{x)pB{y)) - {pAix)){pBiy)), (where p^(x) = ^Tr S{x-A) and similarly pBix)) is, in the 
large N limit, 
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pAB{x,y) 



2y;-2 (.Qg (,Qg 

1 + chu cos{9 + a) 1 — chu cos{9 — a) . ,^ 



' (chw + cos(6' + a)y {chu — cos(6' — a)Y 
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Here Vi(x) = V2(x) = (l/2)/ix^, u = In(^), x = asinO, y = asina. The one point function 
p(x) = {pa{x)) = {pb{x)) is still given by the semicircle law except that the endpoint of its 
support is now a = {^^Y^"^, and (p!0|) is vahd for |y| < a. (|T0|) was computed in ||23 



for P = 2 using a diagrammatic technique and in for /3 = 1,2,4 using the method of 
loop equations 



The local fine-grained correlators in the 2-matrix model for the gaussian ensembles were 
computed in (/? = 1, 2) and |^ {f3 = 4) using the the supersymmetry technique 



For (3 = 2, global fine-grained correlations are known (see also |^) using the method 



of orthogonal polynomials The local correlators near the centre of the eigenvalue 

distribution have been shown to be universal (independent of the form of the potential) in 
T0| (see for numerical evidence and for a physical argument). The smoothed global 



correlator (0) although derived for a quadratic potential reduces, in two different limits, to 
expressions that are universal. The first is when x, y are restricted to the region near the 
centre of the cut; then the expression (see section 5) is the local smoothed correlator 



| 35| which is universal. The second is when /i and c are tuned such that the probability 
distribution has vanishing weight over configurations in which A ^ B, in which case 
(PUP reduces to the universal global 1-matrix result (|^). However, the universality of (|r^) in 
general remains unproven. 

Correlators in the double scaling limit of the (3 = 2 2-matrix model with non-quadratic 



potentials have been studied in [36,37 



3. Quantum Chaos 

The subject of quantum chaos is at present defined as the study of the quantum properties 
of a system that is classically chaotic. Consider a classically chaotic system described by 
some hamiltonian H = H{pi, qi,P2, (I2, •••)• ^ signature of classical chaos is that trajectories 
are extremely sensitive to initial conditions. The question is, what can one say about the 
quantum system? Here we restrict ourselves to questions regarding the spectrum of the 
quantum hamiltonian. In general, the analytic calculation of the energy eigenvalues is a 
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hopeless task. However, there is a substantial amount of data on the discrete eigenvalues Ei 
obtained numerically for model systems and experimentally for real systems. From this data 
one obtains various statistical properties of the spectrum. For example, one can in principle 
determine the correlation functions of the density of eigenvalues {p{Ei)p{E2)...p{En)) , where 
p{E) = J2i S{E—Ei). (( ) is an average over many energy levels of the hamiltonian in question 
keeping the differences Ei — Ej fixed.) In practice one computes various statistical measures 
from the data that are related to these correlation functions (for reviews see 
for example, the distribution of level spacings, (i.e., the probability for the spacing between 
consecutive energy levels to lie in a specified range). The result of these empirical analyses 
is that while the one-point function {p{E)) varies from system to system, the local two and 
higher point correlation functions as well as the level spacing distribution, for sufficiently 
high energy levels not related by any symmetry, exhibit a universal (i.e., system independent) 
behaviour for a wide variety of systems. 

This observed universality of local eigenvalue statistics suggests that there might be a 
probability distribution of energy eigenvalues with universal features characterizing generic 
chaotic systems. It was Wigner's insight in the context of nuclear physics, that the 
statistical spectra of complex nuclei should be obtainable by taking the extreme view that 
one knows essentially nothing about the hamiltonian, that is, the hamiltonian is drawn 
completely at random from the space of all possible hamiltonians. This led him to propose 
gaussian random matrix ensembles for theoretically calculating the eigenvalue statistics, 
wherein the random matrix A of the 1-matrix model in the previous section corresponds to 
a hamiltonian, and correlation functions of the 'observable' p{x) in the matrix ensembles 
are to be compared with empirically obtained correlators of the eigenvalue density p{E) {x 
is identified with E upto a rescaling to be discussed later). Dyson |1| showed that the three 
different ensembles in random matrix theory should correspond to hamiltonians in three 
different universality classes identified by symmetry: (3 = 2 corresponds to hamiltonians 
with no time reversal invariance; if H has time reversal invariance, then (3 = 1, except if the 
system has no rotational invariance and has odd spin, in which case (3 = 4. It was proposed 
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38 Hi, 



that the same statistical hypothesis should apply to complex atoms |^ and to electrons in 
disordered cavities (mesoscopic systems) Later, it was suggested |Q and numerically 
observed that spectra of classically chaotic systems obey level repulsion as in RMT, 
and a proposal was made in |^ that that local correlations in RMT should provide the 
eigenvalue statistics for any hamiltonian system in the strongly chaotic regime. 

These proposals have been tested and verified in an amazing diversity of systems. The 
list includes the empirical spectra of complex nuclei atoms molecules 



and microwave cavities the numerically obtained spectra of chaotic billiards |^ and 



other chaotic systems Q], and analytic calculations from microscopic models of electrons in 
weakly disordered mesoscopic systems . The empirically observed universality mimics the 



universality mentioned for RMT in section 2, in that the eigenvalue density is nonuniversal 
but higher point local correlators and the spacing distribution is universal. For systems in 
which one can tune parameters in the hamiltonian to go from the classically integrable to 
the chaotic regime, agreement with RMT is seen in the strongly chaotic regime, and the 
statistics departs from RMT in the integrable or near integrable domain. 

Recently, the smoothed global correlators have also found applications. An application 
of the universality of is an argument for the universality of conductance fluctuations 



of weakly disordered mesoscopic conductors using the transfer matrix formalism I^T] . Other 
examples include possible application to the spectrum of the Dirac operator in strongly 
coupled QCD and to the quantum Hall effect. 

It seems that while being classically chaotic is usually a good indication of a system being 
'generic' enough for Wigner's statistical hypothesis to apply, there are some exceptions to 
this rule. One class of exceptions are conductors with strong disorder. Classically, these 
systems represent chaotic billiards since an electron can bounce repeatedly from randomly 
placed impurities. In practice their spectrum disagrees with random matrix theory. In this 
case the departure can be explained in terms of a new quantum phenomenon, localization, 
that arises due to strong disorder, which disallows quantum states that are analogues of the 
classical ballistic chaotic trajectories [|^. Other exceptions are discussed in [Q. Therefore, 



until a better understanding of the precise domain of applicability of RMT is reached (see 



for efforts using semiclassical methods |5^), the agreement of energy level statistics 



of classically chaotic systems with random matrix theory remains a fairly ubiquitous, but 
essentially empirical, fact. 

4. The Integrable Model: Static Correlators 

Consider the system of fermions in one dimension with the hamiltonian (the Calogero- 
Sutherland-Moser (CSM) system 



N . . 

H = Eb' + ^"^") + oPiP - 2) E 7 V2- (11) 



The particles are in a harmonic potential and interact with each other via the potential 
with a strength determined by /3 (/3 = 2 corresponds to free fermions). This is an integrable 
model with N conserved quantitities. Its exact ground state wave function is given ||5^ by 



,1 

i<j 

Therefore the expectation value of any position dependent operator F(xi, xtv) in the 
ground state is given by 

(F) = J dxi...dxNi^*Ftp = J dxi...dxNYl\xi- Xjfe-'^^^''^^F{xi,...,XN). (13) 

i<j 

This has the same form as the expectation value of an observable in the Wigner-Dyson 
distribution (§) with V{x) = j^ojx^, = Nq, at /? = 1, 2, and 4, with the position 
coordinates of the particles corresponding to the eigenvalues of the random matrix. Thus 
the matrix model (|l|) describes free fermions for (3 = 2 and interacting fermions for /? = 1, 4. 
In particular, for F we may choose the n-point function of the density operator in the CSM 
model, p{x) = I]j 5(x — Xj). Then, for example the connected spatial equal-time density 
correlator in the 1/r^ model, exactly matches the matrix model and hence the quantum chaos 
eigenvalue density correlator: (p(x)p(?/)) |csm model = N{p{x/ y/N)p{y/ ViV)) Matrix model, with 
the r.h.s. given by (|^) or (||). 
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5. Dynamical Correlators, Perturbed Chaotic Hamiltonians, and the 2-Matrix 
Model 

We now discuss a deeper connection between quantum chaos and the CSM model, established 
through the intermediary of a 2-matrix model. Consider a chaotic hamiltonian H consisting 
of an unperturbed part proportional to Hq and a perturbation proportional to Hi. 

H = Hocosn(f) + Hi^^^. (14) 

is a measure of the strength of the perturbation {Q is a constant, introduced for conve- 
nience, to be specified shortly). Given Hq and Hi, the eigenvalues of H are a function of 
(f). One can think of a 1-d gas of particles whose positions are the eigenvalues and (p as the 
time on which they depend. An interesting fact is that the time evolution of this gas (0 
dependence of the positions of the eigenvalues) is given by a classical hamiltonian which is 
closely related to the CSM hamiltonian with some additional degrees of freedom. This is 
known as the Pechukas gas [Q. We do not discuss this further; instead we focus on the 
'quenched' averages that arise when the spectral statistics of H is discussed using matrix 
model ensembles. In particular we discuss the dependent correlator {pho{Ei)ph{E2))c 
which is of interest in quantum chaos, since it measures the correlations between the eigen- 
values of Ho and if as a function of the strength of the perturbation {ph{E) = J2 ^{E — Ei) 
where E^ are the eigenvalues of H). 

To see how this correlator is connected to the 2-matrix model, note that the prob- 
ability P{H, (j))dH that the full hamiltonian lies in the volume element dH located at 
H (in the space of N x N matrices of the appropriate class given by f3) is given by 
P{H,(f))dH = J dHoPoiHo)PiH,(f)\Ho)dH, where Po{Ho)dHo is the probability that the 
unperturbed hamiltonian lies in the volume element dH^ located at Hq, and P{H, (f)\HQ)dH 
is the conditional probability that full hamiltonian lies in the volume element dH lo- 
cated at H given a fixed unperturbed part Hq. But now P{H,(f)\HQ)dH = Pi{Hi)dHi = 
const. Pi{Hi)dH , where Pi{Hi)dHi is the probability of the perturbation being in a 
volume element dHi located at Hi, and we have assumed that Hi and H come from 
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the same class of ensemble (same whereby, for fixed Hq, dHi = const. dH . Then 
P[H,(f))dH = const. J dHQPo{H())Pi{Hi)dH , where Hi is understood to be given in terms 
of Ho and H by (0). Thus the average of any H dependent quantity will be given by a 
2- matrix integral (the two matrices to be integrated over being Hq and H), in which the 
factor Pi{Hi), when expressed in terms of Hq and H, will involve (0 dependent) terms that 
couple the two matrices together. In particular if we take both probability distributions Pq 
and Pi to be gaussian: Po{Ho) ~ exp(— ^Trifg), Pi{Hi) ~ exp(— |TriJf ), it follows that 

P{H, (j))dH = const. J dHodH exp(-Tr[F(i7o) + V{H) - cHqH]), (15) 

V{x) = l-fxx^, /i = . ^ , and c = (16) 
2 sm \l(p sm iZ0 



It is convenient to rescale Hq = y/NA, H = \'NB, and define the 2-matrix model 
partition function (y) with S = NTt {V{A) + V{B) — cAB). From the symmetry between 
Hq and H in (|T5]) it follows that the two eigenvalue densities are equal, {pa{x)) = {pb{x)) = 
p{x), and are given by the semicircle law for gaussian Pq and Pi. The constant Q is chosen 



to be f2 = \J Ti"^ (3 / 2N which means that /i and c are 0(1) when ~ 0(1). It is evident 
from the above that under the random matrix hypothesis the desired 2-point correlator in 
quantum chaos is given by p2| , |39| , |6l | 



c I quantum chaos 

Npab{E'i,E'^), (17) 

with E[ = N~^Ei. 

Altschuler, Simons, Szafer, and Lee provided evidence based on numerical results for 



chaotic billiards and the Anderson model for the universality of the l.h.s. of (|T^ ||32| , p6 
computed the matrix model local correlator for the r.h.s. of ( p!7D and exhibited agreement 
with numerical results ||26|, and conjectured that the r.h.s. of ( p!7D furthermore equals the 



(iynamica/ density-density correlation function in the CSM model with the perturbation 
strength in the quantum chaos problem playing the role of time in the CSM model (and 
Ei playing the role of the positions Xi of the CSM particles). 
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We now outline the proof of the latter conjecture |^ refering the reader to for 



a detailed review. (See also |^ for related observations and for a proof that holds for 
smoothed correlators only, see [^.) It is not difficult to see that the probability distribu- 
tion P{H,(f)\Ho) = const. exp{—TT[V{HQ) + V{H) — cHqH]) for the matrix elements of H 
which appears in (|T3p is a solution of a Fokker-Planck equation describing the independent 
Brownian motion of the matrix elements of H, where the time t is given by 

log(cos(l]0)) 

t- , (18) 

and the initial condition is that at t = = 0, the distribution is localized on Hq, 
P{H, 0\Hq) = 6{H — Hq). Dyson |^ had shown that when the elements of a matrix exe- 



cute independent Brownian motion, its eigenvalues execute a correlated Brownian motion 
in which they repel each other by a logarithmic interaction, and derived the corresponding 
Fokker-Planck equation for the time dependent probability distribution of the eigenvalues. 
Sutherland showed that the eigenf unctions of this latter Fokker-Planck hamiltonian 
were in one-to-one correspondence with the eigenfunctions of the CSM hamiltonian. It is 
this correspondence, which extends the relationship between the matrix model and the CSM 
model at the ground state level (responsible for static correlators as discussed in the previous 
section) to all the excited states, that is at the heart of the mapping between the dynamical 
correlators of the CSM model and matrix models. One can also map the CSM model to a 



continuous matrix model in which the matrix A is a function of a continuous variable t 64 



For the two-point function the precise mapping is ||6T] 



(p(x, 0)piy, t))c|csMmodei = Npab{x', v'), (19) 

where p{x,t) = J2iS{x — Xj(t)), and Xiit) is the time dependent position operator of the 
CSM model in the Heisenberg picture. The l.h.s. of (p!9|) as a function of x,y,t, and 
parameters u, (3, N is thus expressed as a correlator of the 2-matrix model (|^) with couplings 
p, c expressed in terms of t via eqs. (|T^|T^), u = and the mapping works for the values 



/? = 1, 2, 4 for which the r.h.s. can be defined in terms of the appropriate matrix ensemble. 
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Eq. (|T0|), substituted in ([I9D, therefore gives the smoothed global spacetime dependent 
density-density correlator of the CSM model. If one specializes (|T0|) to the region near the 
centre of the semicircle {x,y <^ a), and for t cj)^ <^ 0{N) one recovers the previously 



known 'translationally invariant' result |35 



(p(x, 0)p(,, t)). ^ -^1^3^^^^- (20) 

This identifies the sound velocity as Vs = vr/?. We mention that for the CSM model at other 
rational values of (3 not covered by known matrix ensembles, local fine-grained dynamical 
correlators have been obtained by the method of Jack Polynomials . See for a proof 
of (H) for all even /3. 

Together with (|17D the equation ([T9|) estabhshes a connection between perturbation 
strength dependent correlators in quantum chaos and dynamical correlators in the CSM 
model via the 2-matrix model. The r.h.s. of ( ]20|) therefore also equals the parametric 
eigenvalue correlator {pHo{^)PH{y))c in the quantum chaos problem H = Hq + (pHi, when 
Hq and Hi belong to the same symmetry class. We emphasize that the mapping from the 
quantum chaos problem to the matrix model is an ansatz that has essentially empirical 
support. On the other hand the mapping between the matrix model and the CSM model 
discussed above is an exact mathematical mapping. 

6. Two Dimensional Quantum Gravity 

Matrix models have long been studied in quantum chromodynamics, where the basic degree 
of freedom is a matrix valued field. The observation that the expansion of SU(A^) 
QCD groups the Feynman diagrams of the theory according to the topological classes of 
2-dimensional surfaces led to a vigorous study of the random matrix models in the large- iV 



limit ||T9|. This has also led to matrix models being proposed as models of 2-dimensional 
quantum gravity and string theory where one needs to sum over all 2-dimensional 

surfaces keeping track of their topology. 

A 2-dimensional surface can be latticized by triangulation, which is a representation of 
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the surface in terms of small, flat, equilateral triangles (all of the same size) joined along 
their edges. In this discrete representation, all information about the local geometry of the 
surface is contained in the local 'coordination number', the number of triangles that meet 
at a vertex. The area of the surface is equal to the total number of triangles (in units of 
the area of each elementary triangle), the length of any boundary is equal to the number of 
edges in that boundary, and the topology of the surface, (given by its Euler characteristic 
X = 2 — 2h — b where h is the number of handles and b the number of boundaries) is 
X = Nq — Ni + N2, where A^o = number of vertices, Ni = number of edges, and N2 = 
number of triangles in the triangulation. 

In a path integral formulation of 2-dimensional quantum gravity one needs to sum over 
all surfaces, assigning a specific weight to each surface. This sum is defined to be the sum 
over all distinct triangulations in the lattice version of the theory |^8[. For example the 



partition function for pure gravity (without matter fields) is defined to be 

00 00 

Z{k, X) = J2 e"^'"""^ E e'^^Z{h, A), (21) 

h=0 A=4 

where Z{h, A) is the number of distinct triangulations of connected oriented surfaces with 
area A, number of handles h, and no boundary, k is referred to as the gravitational constant 
since it multiplies the Einstein action (which is proportional to 2 — 2/i for closed surfaces in 
two dimensions), and A the cosmological constant. 

The above partition function can be obtained from the partition function of a hermitian 
1-matrix model: Z{k,,X) = \ogZ where the Z on the r.h.s. is given by (||), with V{A) = 
(1/2) A'^ — gA'^, N = e'^ and g = e~^. This is apparent upon expanding ([l|) as a power series 
in g, and collecting all the Feynman diagrams which appear with the same power of 1/A^. 
Since every Feynman diagram is a vacuum bubble made of propagators and cubic vertices 
only, its dual diagram (obtained by joining the centres of adjacent loops) is a triangulated 
surface. Since the perturbative evaluation of the partition function involves summing over 
all Feynman diagrams, one automatically gets a sum over triangulations. The weight factor 
appearing with every Feynman diagram containing V vertices, P propagators and L loops 
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is N^~^~^'^g^ , which becomes e''('^''~^i+^2)-AAf2 g^g needed, since every dual diagram has 
Nq = L vertices, Ni = P edges, and N2 = V triangles. The logarithm of Z appears because 
in (^) we need only connected surfaces. The leading large contribution to log Z comes 
from closed connected surfaces with the spherical topology {h = 0, x = '^)- The orthogonal 
and symplectic ensembles lead to sums over non-orientable surfaces (see, e.g., ||69||). 

Consider now the insertion of Tt in the matrix model path integral, i.e., the quantity 
Z(Try4') = / dAe~''^TTA'' , with I a positive integer. Expanding the r.h.s. as a perturbation 
expansion in g, one notes that as before all vertices of the Feynman diagrams are cubic, 
except one (corresponding to the external insertion of TrA') where / propagators meet. 
Thus all the corresponding dual diagrams are triangulated surfaces as before, but with one 
elementary plaquette which has / sides. This distinguished plaquette can be considered to 
be absent, i.e., the triangulated surface can be thought of as having one boundary of length 
/. Thus the insertion of TrA' in the path integral has the interpretation of inserting a 'loop' 
of length / on the surface ||70|. Similarly the insertion of TrA'TrA™ can be seen to produce 



two loops of size / and m on the surface. W{z) = {l/N)J2^=oi'^/^"'^^)'^''^^"' is called the 
loop operator because it is the generating function of loops of all lengths; z is the fugacity 
for the length of the loop. Thus correlation functions of the type 

{W{z,)...W{zn)) = ^JdA e-''W{zi)...W{zn) (22) 

are of direct physical significance in 2-d gravity and string theory since they determine sums 
over surfaces with n boundaries of arbitrary length. 

Now these correlation functions are closely related to the eigenvalue correlations in the 
matrix model because of the simple identity p{x) = {i/27c) lim^_,Q[W {x + ie) — W{x — ie)]. 
Thus for example p{x) is given by the discontinuity across the cut in W{z) and the two 
point function by 

(p(x)p(y)) = \im,^o[{W{x + ze)W{y + ze)) + {W{x - ze)W{y - ze)) 

-{W{x + ie)W{y - it)) - {W{x - it)W{y + it))]. (23) 
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Therefore, in view of the correspondence discussed in the previous sections, a single math- 
ematical quantity in the 1-matrix model (the r.h.s. of (|2^)) simultaneously gives three 
different things: correlation functions of loop operators in 2-d gravity, eigenvalue correlators 
in quantum chaos, and static spatial correlators in the CSM model. 

The 2-matrix model defined by with a cubic potential V{A) = ^A^ — gA^ maps onto 
2-dimensional gravity coupled to the Ising model ||7l| . The two matrices A and B correspond 



to the two spin configurations of the Ising model and the coupling constant c in the action 
S{A, B) provides the extra Boltzmann factor when adjacent spins are unequal. 

7. Loop Equations and Eigenvalue Correlators 

Loop equations |7^Jl9[| are identities relating different correlation functions of the loop oper- 
ators. There is a hierarchy of identities through which successively higher point correlation 
functions get coupled to the lower (e.g., one or two point) correlation functions. In the large 
N limit, the connected higher point correlation functions are associated with higher powers 
of as compared to the disconnected parts; this causes the infinite hierarchy of loop 
equations to break up into finite subsets of algebraic equations, from which all correlation 
functions can be solved for by an inductive procedure. 

As a simple example of this procedure in the 1-matrix model consider the identity 

where A is a real symmetric matrix, S = NTiV{A), V{A) = (1/2) fxA"^ (gaussian orthogonal 
ensemble, /? = 1), and i and j are not summed over. Expanding the action of the derivative 
on its argument and summing over i and j one gets the identity (for details see [^ ) 

= -fiizWiz) - 1) + ^{Wiz)Wiz)) - ^§-Wiz). (25) 

This is an exact loop equation where the one and two point correlators are both present. Now 
decompose {W{z)W{z)) into its disconnected and connected parts: (W^zY) = {W{z)y + 
{W{z)W{z))c. In the large N limit {W{z)) = W{z) is 0(1) and {W{z)W{z))^ is O(^) (see 
below). Then to leading order in 1/A^ we can suppress the latter as well as the last term in 
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(p5|), obtaining the closed large loop equation = —fi{zW{z) — 1) + {l/2)W{zy with the 
solution W{z) = {2/a^)[z — ^ z^ — dP\ which has a branch cut between z = ±a, = 2(3/ jj,. 
Using the identity p{x) = {i/2TT) lim^^Q[W{x + ie) — W{x — ie)] in the above solution for 



ically, E ~ 0{N'^) (e.g., for the gaussian case when the endpoint a = y2/?//i is 0(1), all 

the eigenvalues Aj are in an 0(1) region around the origin, and then both the terms in 
(I) are 0{N^)). Thus to leading order F = InZ is 0{N'^), say F ^ N'^T, where T is a 
function of Qn only, independent of A^. From this it follows that {TtA"'TtA"^)c = dg^dg^T is 
0(1), while the disconnected part (TrA")(TrA"^) is 0{N^) since (TrA") = -Ndg„T is 0{N). 
Alternatively one can understand factorization geometrically. As indicated earlier, the 
dependence coming from various types of surfaces is proportional to A^^, where x is the 
Euler characteristic of the surface. Let us restrict our consideration to surfaces with h = 
(no handles) which will always give the leading large A^ contributions. {W{z)W{z))c, since 
it contains two insertions of the loop operator, corresponds to connected surfaces with two 
boundaries whose Euler characteristic is therefore x = 2 — 2/i — 6 = 0, while the disconnected 
part of {W{zy), namely {W{z))^, corresponds to a pair of surfaces each with one boundary 
(and hence each with x = 1), thereby representing a total x = 2. Thus the disconnected 
part is a factor A^^ larger than the connected part. 

To obtain the 2-point function of the loop operators, start from the identity 



Expanding the action of the derivative on its argument yields an identity connecting 2 
and 3-point functions of the loop operator. Expanding that identity into connected and 
disconnected parts as before, one finds that terms proportional to ( ^5]) exactly cancel, while 
the connected 3-point correlator piece is suppressed by factors of 1/A^ or 1/A^^ compared to 
other terms. In the large A^ limit one gets the closed loop equation 



W{z), one immediately reproduces the Wigner semicircle law (|^). 

To prove 'large A^ factorization' one observes that in (^) when Qn are 0(1), then gener- 





(26) 



= {W{z) - fiz){W{z)W{w)), + 



1 d W{w) - W{z) 



(27) 



A^^ dw w — z 
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which yields upon simphfication 

^ ^ N^2P{z-wy{l-3^W^iz)){l-S^W^iw)) ^ ' 

We have outhned the proof of this equation for the gaussian orthogonal ensemble, = 1. The 
same steps repeated for the other two ensembles yields the (3 dependence displayed above. 
The difference enters at the point where action of the derivative is expanded. Concretely, for 
(3=1, 1^ = 6ikSji + 6ii6jk whereas for (3 = 2, |^ = 6ikSji, and a slightly more complicated 
expression exists for /3 = 4. The difference in these expressions simply reflects the different 
number of independent degrees of freedom in the three ensembles. The extension of the 



method of loop equations to the /5 = 1,4 ensembles was discussed in From (^) and 
(p3|) it immediately follows that the 2-point connected correlator of the eigenvalue density 
operator is given by (P). In ref. |jl3[ (pH]) is derived for an arbitrary potential from the 
method of loop equations for (] = 2. That also establishes the universality of (P) for the 
unitary ensemble. We expect that this proof of universahty of ( pSj) would also extend to the 
other ensembles. 

For the 2-matrix model, since the two matrices are coupled, a larger number of identities 
are needed to get closed equations. A procedure similar to the one outlined above yields for 



the three gaussian ensembles the global correlator 

1 rn"^ 1 

{WAiz)WBiw)), 



2fi(3{l-^W{z)W{w))^ 
, W{zf W{wf ^ ^29) 



with a = (4^)^/^. (0) follows from (|9|) upon using (|23|). This method bypasses the 
problem of angular integrations. For other applications of loop equations in hermitian, 
non-gaussian 2-matrix models, see |37[] . 

Finally we remark that (^), ® , (^Of), and ( p!0D provide 'smoothed' correlation functions 
because we first compute loop operator correlation functions for z and w far from the cut 
which automatically averages over the oscillations on the scale of the eigenvalue spacing. 
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8. Summary 

We have described some elementary mathematical properties of random matrix models and 
their connection with three distinct areas: two-dimensional quantum gravity, the in- 
tegrable model, and quantum chaos (the last connection includes applications to nuclear 
physics and mesoscopic systems). The loop operators are a link between these areas since 
they provide eigenvalue correlations in quantum chaos and the CSM model on the one 
hand, and also have a direct geometric significance in quantum gravity and string theory. 
The method of loop equations is a powerful method for calculating smoothed global cor- 
relators. This method treats all ensembles on the same footing because the procedure for 
deriving the loop equations remains the same. The only distinction between the ensembles 
enters at the level of counting the degrees of freedom, that is, in the expression for 
The extension of loop equations to the orthogonal and symplectic ensembles should also find 
application to sums over unoriented random surfaces. 
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